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Hypernormal distribution theory: Analysis of the set of extreme random variables 
models 


Abstract: The analysis of the set of extreme random variables models is still an extremely topical topic in 
many areas of mathematical research in the theory and practice of managing production processes due 
to its specificity and great interest in finding an expectation and stability indicators set studied in practical 
economics. Calculations of applied mathematics help to determine tentatively possible boundary 
parameters of various models, t.c., expectations, despite the fact that theoretical calculations do not have 
a direct association with practical data. Nevertheless, the consideration of extreme models of extreme 
random variables is still relevant in many areas of science and industry. The study subject was the 
hypernormal distribution theory. The study object was a set of extreme random variables models. The 
study purpose was a comprehensive analysis of many models of extreme random variables. To achieve 
the purpose and solve the tasks formulated on its basis, empirical, analytical and comparative methods of 
data analysis and the method of mathematical modelling, which contributed to the study of the materials 
presented in this article, were used. In the study course, materials from the works of such leading world 
experts in extreme value theory and programming K. Beck, M. Fowler, L. Tippett, E. Gumbel, K. Auer, 
R. Miller, and Scott W. Ambler and researchers as V.L. Khatskevich, B.V. Gnedenko, V.A. Akimov, V.A. 
Bykov, E.Yu. Shchetinin, K.M. Nazarenko, L.P. Kvashko, A.S. Losev, V.S. Mikhailov, V.A. Popov, E.R. 
Smolyakov. In the study course, the definition of an extreme value within the framework of the theory 
was refined, the typology of the distribution of maximum values was analysed, seven theories of the 
hypernormal distribution were identified and their proofs were presented, and practical examples of the 
application of each theory were given. The practical significance of the study of extreme random variables 
models in various areas of industrial human activity was confirmed. The materials of the study can be 
used in the widest range: from application in risk management of industrial production to predicting the 
probabilities of natural phenomena, which makes it possible to prevent significant economic and social 
losses of society, as well as make a tangible contribution to programming the probabilities of the 
development of the society of the future. 
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Teopua rumepHopMasbHOro pacnipeAecacHHA: AHAAH3 MHO2KECTBA MOACACH 
SKCTPeMaAbHBIX CAyUAMHbIX BEAMUHH 


Annomayua: AHaAW3 MHO?KECTBA MOACACH SKCTPEMAABHBIX CAYYAMHIX BEAMUMH AO CHX Op ABAAeTCA 
Kpalive akKTyaAbHOH TeEMOIt BO MHOrMX OOAACTAX MaTeMaTH4eCKHX HCCACAOBaHHI TeOpUu MW WpakTHKH 
ylipaBACHHA IIPOH3BOACTBCHHBIMH lIpomeccaMu B CHAY cBoel cnenmcbu4HocTu uu OOABLIOTO uHTepeca 
K TIOHCKy MHO?KCCTBa TlOKa3aTeAecit O?/KMUAACMOCTH CTaOMAbHOCTH, HMICCACAYVCMBIC B rpakTHueckou 
SIKOHOMHKE. Pacuérsr TIPpHKAaAHOHW MaTeMaTHKH TIOMOlartoT OIIPCACAATb OpHCHTHPOBOUTHO BO3MO?KHBIC 
MOrpaHHatHple TlapaMeTpbl BCCBO3MO?KHDIX MOA€AeH, T.¢., OowKMAAHUH, HECMOTPA Ha TO UTO 
TeOpeTuaeckKHe pacuérer He UMCFOT rIpaMon accoluHalun C TIpakTHueCKHMMU A@aHHbBIMH. Tem He MecHee, 
paccmMoTpeHue IKCTPCMaAbHbIX MOACACI IKCTPCMaAbHbIX CAY4aHplx BCAHUHH AO CHX mlop ABAACTCA 
aKTYAABHDIM BO MHOITMWX OOAaCTAX HaykKu uu TIPOMBIMIACHHOCTH. IIpeametom MCCACAOBaHHA ABAHAACb 
TeopemMa THIeC¢pHOPMaAbHOTO PpaciipcAcAcHHa. OOpextTom MCCACAOBaHHA ABAAAOCB MHO?KECTBO MOACAcIt 
SKCTPCMaAbHBIX CAyU4aHHBErx BeCAMUHH. Lleanro AaHHOTO HMCCACAOBaHHA ABAHACH KOMIIACKCHBIM aHaAH3 
MHO?KCCTBA MOACACK IKCTPCMaAbHBIX CAyUdaHHBErx BeAMUHH. AAS AOCTVWKCHMA TIOCTaBACHHOM WeAu UW 
peltennA ccbopMyAHpOBaHHBIX Ha eé OCHOBaHMM saAat MCIOAB3OBaAMCb IMIMpHueckHH, 
aHaAHTH4ecKHM iu CpaBHHTeAbHEIM MCTOADI aHaAH3a AaHHBIx uu MCTOA MaTeMaTU4eCKOrO 
MOACAMPOBAHHA, KOTOPbIe CIOCOOCTBOBaAH HCCACAOBaHHIO MaTepHaAOB, UPeCACTaBACHHBIX B AaHHOIt 
cTaTbe. B XOAC HMCCACAOBAHHMA OBIA MCIIOAB30OBaHBl MaTepHaAbl TPYAOB TaKHUX OTCUCCTBCHHbIX 
uccAeaopateAen kak B.A. Xarnkesuy, b.B. PHeaenxo, B.A. Axumos, B.A. Barxos, E.fO. Werunun, K.M. 
Hasapeuxo, A.J]. Ksanrko, A.C. Aoces, B.C. Muxaiiaos, B.A. [lonos, 9.P. Cmoapaxos, a TamKe 
MaTepHaAbl BeAYUIHxX 3apyOeKHBIX CIeHWHaAWMCTOB B oOAacrH TeOpHu IKCTPCMaAbHbIxX 3HaudcHM iu 
mporpamMuposanua K. Bexa, M. Mayaepa, A. Tunnerra, 0. Tym6eas, K. Ayapa, P. Muaaep, u Cxorra 
B. DmOaepa. B xoAe uccAcAOBaHHA OBIAO YTOUHEHO ONPeACACHHE SKCTPEMAABHOM BEAWYHHBI B paMKaX 
TCOpHH, WpoaHaAv3upoBaHa THITOAOTMEO PpactiIpeAcAcHHA MaKCHMaAbHbIxX BCAMUHH, OIIPCACACHBEIL Ce€Mb 
TeOpeM THMeCpHOPpMaAbHOTO pacipcAcAcHuA uw TIpCACTaBACHBI UX AOKa3aTCAbCTBa, a TakKoKe AaHbIi 
UpakTH4eckHe TIPHMEPppbI UpHMcHeHHA KaoxKAOM M3 TeOpeM. Tem CaMBIM, ObrAa TIOATBep7#KACHa 
UpakTHu4eckaA 3SHaAYUMMOCTB HMCCACAOBaHHA BapHaTHBOB IKCTPeCMaAbHBIX MOACAcII IKCTPCMaAbHbIX 
CAyudaHHBEIXx BCAHMUHH B pas AM4HbIX OOAaCTAX MHAYCTpHaAbHon ACHTCABHOCTH AMCAOBCKA. Matepraapr 
AaHHOTO HCCACAOBaHHA MOTYT OBIT MCIIOAB3OBaHbI B CaMOM INMPpOKOM CHeKTpe: OT TIPHMCHCHHA B 
oOaAacTH PHCK-MeCHeCAKMCHTa TIPOMBIITACHHOTO IIPOH3BOACTBa AO TIpeAcka3aHHA BepOATHOCTEHM 
TIpMpOAHBbIx ABACHHH, 4TO IO3BOAAHCT TIpecAylpcAHTb SHaAYUHMTCABHBIC SKOHOMUMYCCKHE MW COMUMAABHDHIC 
TWloTepH oOmectBa, a TaKoKe BHCCTH OLIYTHMbIM BKAaA B IIpOrpaMMUpoBaHHe BepOATHOCTeH pa3sBuTuAa 
oOulecrBa OyAylero. 


Kutouesoie cno6a: IKCTPCMaAbHBIC CAYWaHHBIe BCAMUMHDI, THIe€PpHOPMaAbHOoe paciripeAcaAcHue, TCOpHA 
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Introduction 


SKCTPCMaAbHBIX 3HauecHHH, SHTPONHA. 


The analysis of the set of extreme random variables models is still an extremely topical topic 
in many areas of mathematical research in the theory and practice of managing production 
processes due to its specificity and great interest in finding an expectation and stability indicators 
set studied in practical economics. Calculations of applied mathematics help to determine 
tentatively possible boundary parameters of various models, i.e., expectations, despite the fact 


that theoretical calculations do not have a direct association with practical data. Nevertheless, 


2 


the consideration of extreme models of extreme random variables is still relevant in many areas 
of science and industry. 

The study subject was the hypernormal distribution theory. 

The study object was a set of extreme random variables models. 

The study purpose was a comprehensive analysis of many models of extreme random 
variables. 

Based on the study purpose, the following tasks were formed: 

— clarify the definition of an extreme value within the framework of the theory; 

— analyse the typology of the distribution of maximum values; 

— define hypernormal distribution theories and present their proofs; 

— give a conclusion on the practical application of the evidence base of the hypernormal 
distribution theories. 

To achieve the purpose and solve the tasks formulated on its basis, empirical, analytical and 
compatative methods of data analysis and the method of mathematical modelling, which 
contributed to the study of the materials presented in this article, were used. 

In the study course, materials from the works of such leading world experts in extreme 
value theory and programming K. Beck, M. Fowler (Beck ¢> Fowler, 2001), L. Tippett (I7ppett, 
2013), E. Gumbel (Gumbel, 2012), K. Auer, R. Miller (Auer ¢ Miller, 2007), and Scott W. Ambler 
and researchers as V.L. Khatskevich (Khatskevich, 2013; Khatskevich, 2020a;,; Khatskevich, 2020D), 
B.V. Gnedenko (Guedenko, 1943), V.A. Akimov, V.A. Bykov (Akimov et al, 2009), E. Yu. 
Shchetinin (4Aimov et al, 2009; Shchetinin & Nazarenko, 2008), K.M. Nazarenko (Shchetinin ¢ 
Nazarenko, 2008), L.P. Kvashko, A.S. Losev (Kvashko ¢ Losev, 2013), V.S. Mikhailov Mikhailov, 
2012), V.A. Popov (Popov, 2073), E.R. Smolyakov (Swolyakov, 2077). 


Materials and methods of research 


Auxiliary information and basic definitions 

The theory of extreme values is a branch of the science of statistics, which aims to study 
extreme deviations from the median of probability distributions, ie., an assessment of 
phenomena based on an ordered selection of probability parameters for the most extreme events 
ot processes. The concept of extreme value theory was introduced by Leonard Tippett (17ppez, 
2013) in the first quarter of the 20" century and became the basis of many studies that have been 
going on for about 100 years. At that time, his research was based at the British Cotton Research 
Association, where he worked on strengthening the cotton thread. In his research, L. Tippett 
postulated that the strength of a thread is determined by the strength of its weakest fibers. He 
obtained three asymptotic limits that clearly described the distributions of extrema that 
considered independent variables (T7ppett, 2072). It was the study that became the starting point 
in applying a qualitatively new approach to calculating extremeness in production and economic 
indicators. In the future, ED. Gumbel codified this theory in his work Statistics of Extremes 
(Gumbel, 2012). There he gave the distribution concept, which now bears his name. In the second 
half of the 20" century, the results obtained were significantly expanded and began to consider 


insignificant correlations between variables. Strong correlations of the order of dispersion began 


to be actively studied already at the beginning of the 21“ century with the use of artificial 
intelligence and a neural network. 

With the potentially high probability of extreme manifestations, there naturally becomes an 
increased risk of redundant programming as a form of agile software development 
methodologies. The authors of this methodology are such prominent scientists today as Kent 
Beck, Ward Cunningham, Martin Fowler (Beck ¢» Fowler, 2007) and others. Kent Beck pioneered 
the development of the methodology for the Chrysler Comprehensive Compensation System 
project (Beck, 2003). The goal was to apply theoretical methods and develop new and modern 
software for those times. As a result of the development, it was possible to raise and develop 
technology and programming at a new qualitative level. It should also note that it was in extreme 
programming that a departure from the long-term process of creating programmes was 
determined, which consisted in the fact that instead of one-time planning, analysis and design of 
a system for the calculated course of events, specialists now implement these operations in a 
phased complex during development. 

The analysis of extreme values plays an important role in the study of many phenomena and 
in solving applied problems of the complex systems reliability and efficiency, structural 
mechanics, the theory of stability, dynamic strength, etc. Consideration of absolute extrema will 
begin with consideration of the maximum: 

US Max Ky, Xap xcig Xn) 

Values from a set of n random variables (random sequence). If all components of the sample 
X1,..Xp, ate independent and equally distributed random variables, then the distribution 
function of the largest F, (X) value is determined as follows: 

F(X) = P{U < x} = P{X, < X,X2 < X,...,X_ < X} = F™(X) 
where F(X) is distribution function of the original random variable. 

If V is the minimum value of a random variable from a set of 7 random variables: 

V SMa A, 2 Xe): 

And if the components of the sample X,,...,X, are independent and equally distributed 
random variables, then the distribution function of the smallest Q,(X) value is determined 
similarly: 

0,(X) = PIV < xf=1— P{X, = X,X, =X, ...,X, =X} =1—-—[1—F"@)]". 

Thus, extreme values distributions can be derived from the exact original distribution. In 
reality, the analytic properties of the original distribution are rarely known. This leads to the need 
to use the principle of maximum distribution and determine on this basis the extreme 
distribution of the extreme value (maximum or minimum). From a mathematical viewpoint, the 
maximum principle’s application uncertainty leads to the solution of extremal (variational) 
problems with organic ones, determined by the form of setting the probabilistic characteristics 
and the range of random variable values. 

According to the distribution, which has the greatest entropy under certain restrictions, is 
called extreme. Next, a brief description of 8 types of extreme distributions of extreme random 
variables with a certain degree of universality will be given. The common thing in the formation 
of such models is the definition of the Euler-Lagrange equations of variational problems, 


considering the specifics of specifying information about the initial random variable and allowing 


meaningful interpretation. The most obvious and, most importantly, the most practical 
application-oriented is the statistical interpretation of extremal distributions in terms of the 
theory of order statistics, the subject of which is the study of the properties and applications of 
otdered random variables and functions of them. To this end, we present some auxiliary 
information from the theory of order statistics. The source material for statistical analysis, 
obtained as a result of a simple random selection from the general population, determined by 
the random variable X, is a sample of a finite size n: 
Nap ABA 

A sequence of sample values ordered by magnitude x? < X. Me < x”) is called a variation 
series. 

If the initial distribution of the general population is characterised by the mathematical 
expectation m and the variance 6”, the distribution of the rightmost member of the variation 


series is 
Pix < x} = F(X). 

It should let agree to call the distribution that delivers the entropy maximum an extremal 
distribution of type 1. An extremal distribution of type 2 is a limiting n —> 00 distribution of type 
1. If the initial distribution of the general population is characterised by only one mathematical 
expectation M, the distribution of the rightmost member of the variation series is 

Pix < x} = G,(X). 

It should let agree to call the distribution that delivers the entropy maximum an extremal 
distribution of type 3. For the case n > 00 the distribution function G,(X) degenerates into a 
type 4 distribution function G, (X). 

In a similar way, we introduce into consideration the extremal distributions of the minimum 
values. If the initial distribution of the general population is characterized by the mathematical 
expectation m and the variance 6”, the distribution of the rightmost member of the variation 


series (the minimum of the random sequence): 


Pix < x} = Q,,(X). 
It should let agree to call the distribution that delivers the entropy maximum an extremal 
distribution of type 5. 
Asn — © it should say that an extremal distribution of type 5 degenerates into an extremal 
distribution of type 6. If it assumes that the original random variable is characterised only by the 


mathematical expectation m, the distribution of the leftmost member of the variational series 
Pix < x} = R,,(X), which provides the maximum entropy, it should agree to call the 


extremal distribution of type 6. For the limiting case n > ©0 it is useful to introduce a type 8 
distribution. Thus, extreme distributions of types 1-4 are distributions of maximum values, and 
distributions of types 5-8 are distributions of minimum values of random sequences, extreme 
distributions of types 1, 2, 5, and 6 are distributions of extremes of random sequences of 
independent and identically distributed random variables of the general population, the 
distribution function of which is unknown and is characterised only by mathematical expectation 


and variance, extreme distributions of 3, 4, 6, and 8 types are distributions of extrema of random 


sequences of independent and identically distributed random variables of the general population, 
characterised by only one mean, extreme distributions of 3, 4, 7, and 8 types are asymptotic 
(limiting) distributions. In a compact form, the main notation and definitions are presented in 


the appendix (Tabie 7). 


The principle of maximum values distribution 
An extreme distribution of type 1 (hypernormal distribution) is a continuous distribution, 


the probability density of which is the solution of the differential equation: 
n-1 


no*[Fy(X)] » F(X) + &X — m)F, (X) = 0] () 
where m and o* are the mathematical expectation and variance of the set of initial random 
variables. The nonlinear differential equation (1) satisfies the natural boundary conditions 
F,(—c0) = 0, F,, (00) = 1 (2) and is completely determined by the first two moments (m and 
0”) of the original random population and the sample size n. 

The hypernormal distribution corresponds to the distribution function F (X) of the original 
random variable, determined by solving the following differential equation with the same 
boundary conditions 

no? FE, (X)[F(X)]""102n(n — 1) F2(X)[F(X)]"-* + (x — m)E,(X) = 0 (3) 
F(—0o) = 0, F(o) = 1 (4) 

In the appendix, the figures 1-10 (Figure 1; Figure 2; Figure 3; Figure 4; Figure 5; Figure 6; Figure 
7; Figure 8; Figure 9; Figure 10) show graphs of functions and numerical characteristics (expectation 
and variance) of the hypernormal distribution for integer parameters n from 1 to 10, obtained 
as a result of solving the nonlinear boundary value problem (1), (2). The calculation of the 
functions F,(X) is made for standard conditions (for the scale parameter o = 1 the shift 
parameter m = 0). 

For large values of the argument (x >m+30) the hypernormal distribution 
asymptotically approaches the normal distribution with density 


1 
f@)= Jizoun’ ©) 
With an extremal distribution of type 2 (hypernormal distribution), the random variable 


poe has a limiting (for n > 10) hypernormal distribution if its quantile function has the form 


Xp =M+ ovV2n,/—E; (In P) (6) 
where E;(In P) is integral exponential function, whose argument is the natural logarithm of the 
probability p = F(x). 
The second table 2 (Tab 2) presents the values of the function of the limiting hypernormal 


distribution F(x), whose argument is the value: 
x-m 
ovn- 


The mathematical expectation and variance of a random variable is determined by the 


x= 


formulas: 


E|x@| = m+ ov2n J) /—E, Un P)dp = m + 0.69676V2n, 


2 
D = 2no? J [-E:in P)dp — | f. En P)| dp| = 2no(2 — 0.6067?) = 3.0292n0”. 


An extremal distribution of type 3 is a continuous distribution whose probability density is 
a solution to the differential equation: 
mnlGy (x)1G, (x) + G(x) = 0 (7) 
where m is mathematical expectation of the set of initial random variables. 
Differential equation (7) satisfies the boundary conditions 
Gy(0) = 0, Gy(c0) = 1 (8) 
The density is determined by the parameters m and n (sample size). 
An extremal distribution of type 3 corresponds to a quantile function that displays 
DP = On (x)bx 
Xp = mp2F,(L nn 1; Vp) (9) 
where F,(A; B; y;Z) are Gaussian hypergeometric function. 
Using the representation of the Gaussian hypergeometric function as a series, it is 


convenient to represent the calculation formula for function (9) in the form: 
ve) ane 
x = mp Lr-o7P” (10) 


The mathematical expectation of a random variable xO is determined by the dependence: 


(m)] 00 A 
E|Xn | = mn Yr=o (n+r)(2ntr) (11) 
The value of the sum S), series (Table 3): 


Sn = OGG 

An extreme distribution of type 4 is a continuous limiting (n > 00) extreme distribution of 
type 3. The probability density of this distribution G,,(x) is the solution of the differential 
equation: 

MNGoo(X)Goo(X) + Goo(x) = 0 (12) 

with boundary conditions G,.(0) = 0; G,.(00) = 1. 

An extreme distribution of type 4 corresponds to a quantile function that displays: 

P=G(x)Bx, 


X=- art E(In P) (13) 


The mathematical expectation of a random variable x is determined by the dependence: 


E|X,” | = ax [= Jp Ei(in P) dp| = (14) 


ml1n2 


n+1 
In— 
n 


Distribution of minimum values 
An extremal distribution of type 5 is considered to be a continuous distribution whose 


probability density is a solution to the differential equation: 
n-1 _ 
no*[1 — Q,(x)I™ Q,(x) + (x — m)Q(x) = 0 (15) 
where Mm and o7 are mathematical expectation and variance of the set of initial random variables. 


Nonlinear differential equation (15) satisfies the natural boundary conditions: 


Q,(-%) = 0, Q, (0) =1 (16) 


and is completely determined by the first two moments (m and o%) of the initial random 
population and the sample size n. 

The extremal distribution of type 5 corresponds to the distribution function F(x) of the 
original random variable, determined as a result of solving differential equation (15) with 
boundary conditions (15). 

In the appendix, the figures 11-20 (Figure 11; Figure 12; Figure 13; Figure 14; Figure 15; Figure 
16; Figure 17; Figure 18; Figure 19; Figure 20) show graphs of functions and numerical 
characteristics (mathematical expectation and variance) of an extremal distribution of type 5 for 
integer parameters n from 1 to 10, obtained as a result of solving a nonlinear boundary value 
problem (15), (16). Calculation of the functions Q,(X) is made for standard conditions (for the 
scale parameter 0 = 1 the shift parameter m = 0). 

For the values of the argument x <m-— 30, the extreme distribution of type 5 
asymptotically approaches the normal distribution with density: 

fe) =p Re Soa 

As n > © the extremal distribution of type V approaches asymptotically the extremal 

distribution of type 6, whose quantile function has the form: 
X, =m —ov2n J—E(n(1 — P)) (17) 
An extreme type 7 distribution defines a continuous distribution, or density, whose 


probabilities are the solution of the differential equation: 
n-1 oe 
mn[1 — R,(x)]™ R, (x) + RF, (x) = 0 (18) 
where m is mathematical creation of a set of initial random variables. Differential equation (18) 
satisfies the boundary conditions and is completely determined by the first parameters m and n: 
R(0) = 0,R, (=) = 1419) 
An extreme distribution of type 7 corresponds to a quantile function that displays: 
p = R,,(x) Bx; 
mn dom 
Xp = [1-(—p)]* 20) 
(n) 


The mathematical expectation of a random variable X,, 


EX] == a1) 


is determined by the dependence: 


As n > ©, the type 7 extremal distribution ee approaches the type 8 extremal 
distribution. 
An extremal distribution of type 8 is considered to be a uniform distribution with a 
distribution function: 
xm tx<m 
And mathematical expectation: 


E|x@| = © 23) 


Study results and discussion 


This section contains a presentation of the most significant results of the study of the theory 
of extremal distributions of extremal random variables. 
The differential equation (I) defining the hypernormal distribution function F,(%) is the 
Euler-Lagrange equation of the following variational problem: 
a — fr fn(x) In(x)dx > max (24) 
{o., fn@)dx = 1 (25) 
J". f dx = 1 (26) 
(af ex= 007) 
fr — m)? f (x)dx = 0? (28) 
F, (x) = F™(x). (29) 
The validity of this assertion follows from the proof of theories 1 and 2 below. 
Theory 1. Let X is a random variable with density f(x) > 0,x € (—00, 00), F,(X) is the 
distribution function of the extreme member of the variational series constructed from a sample 
of a finite size n from the general population defined by the random variable X. Let, further, the 


first two central moments of the random variable X: 


i [_sfeoar 


= if (x — m)?f (x)dx. 


Then the entropy maximum is reached on a distribution that satisfies the differential 
equation (1). 

Proof. To do this, it is necessary to find the function f(x) and Fn(x), that ensure the 
maximum of functional (24) in the presence of holonomic constraint (29) and under 
isoperimetric conditions (25)-(28). According to the well-known theories of the calculus of 
variations, the Lagrange multipliers A(x),U9,U,,Uz is introduced and the Euler-Lagrange 


equations for the extended function is composed. So, if: 


a —Infn(x) - 1, 


o = -A(x)n[F(x)]"1, 


OL 
af =a A(x), 
the Euler-Lagrange equation for the extended function has the form: 


d Infn(x) _ 
ae + A(x) = 0, 


A(x)n[F (x)]"-1* + vy + 2u2(x — m) = 0. (30) 
The last equation, taking into account (30), can be written as: 


n PLO p(x) = F(x) [vy + 2v2(x — m)]. 


After substituting F(x) = [F,(x)]", the extremal equation in the considered variational 


problem has the form: 
n-1 


n[F,(x)] = F,(x)[v, + 2v2(x — m)]. 31) 


or 


nk, (x) [Frey] + n(n — TF)" 7 F2 (x) = FO) [vy + 202(x — m)]. 32) 

It should integrate equation (32) over the domain of the distribution function F (x), applying 
integration by parts to the first term. Due to certain properties of the distribution function and 
boundary conditions, one can make sure that the integral of the left side of equation (32) will be 
equal to 0, and the Lagrange multipliers v, and U2 will be related to the mathematical expectation 
by the following final relation: 

UV, + 2u2 — 2v2m = 0. 

It follows from here that uv; = 0. Multiplying the left and right sides of equation (32) by the 
independent variable and integrating the resulting equation in a similar way, one can find the 
second final relation connecting the factor with the mathematical expectation k variance. 

Indeed, since: 

SOF?) FO) + (n — 1) FP? (x) F2 (x) dx" (x) F (x), 
then integrating the left side of the new differential equation obtained as a result of multiplication 
by the independent variable gives the following result: 
SO [Fe-1 Oo) Fx) + (n — 1) F"-?2 (x) F2 (x) dx = xF (x) F (x), 
00 


(LF ODEO) + a - DE") F20)] dx = xFOOF@) LS - 


—0Oo 


PGF = —2 


(the first term after the disclosure of uncertainty gives 0). 


Thus, 2u,(m? — o7)2v,m? = —1. Hence it follows that vz = Fe 

Substitution of the Lagrange multipliers vu, and v2 into differential equation (31) makes it 
possible to verify the validity of differential equations (1) and (3). The theory has been proven. 

Many problems of evaluating the efficiency of complex systems and probabilistic analysis 
of complex processes can be formulated in terms of the theory of order statistics and are related 
to the study of extreme values. As an illustration of the foregoing, two examples are given below 
that require such an approach. 

Example 1. The average time to prepare a product for use is 100 minutes, the standard 
deviation is 10 minutes. Five departments simultaneously began preparing for the shipment of a 
batch of five products. Find the probability of preparing the entire batch of products by the time 
T = 120 min. 

Solution. Let F(t) is the distribution function of the product preparation time for use (note 
that in the problem statement the product preparation time distribution law is not specified). Let 
t, is the time of product preparation by the first subdivision, and, u tz is by the second one, and 
so on. 


Then the sequence f(y, tz, ..., ts is a sample of independent identically distributed random 
variables. If to arrange this sequence by the value i < t©) ee to”, the distribution of 


the extreme member of the variational series i) Fey = Fry = P{t® < T} determines the 


probability that the random variable i will be less than the number T. 
According to the condition of the problem, it is required to find 


F,(t)P{t© < 120}. 
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In this case, it is worth noting that such a schematization of a probabilistic experiment 


suggests using the hypernormal distribution function with parameters n = 5, to calculate the 
probability P{t < 120} 


Passing to normalised standard values, it is obtained: 


_ Tom _ 120-100 _ 59 
oO 10 
According to the graph for F5(x) at x = 2.0, it is found: 
F,(2.0) = 0.75. 


Example 2. “Model of collection on alarm”. There are eight individual means (units) of the 
same type (ships, aircraft, etc.) located at the initial moment at the base (airfield, etc.), which, 
according to the alarm signal, should arrive in a given area. The average time of moving one 
object to this region is m = 3 hrs, the standard one is 0.5 hr. In how many hours it is necessary 
to give an alarm signal so that all eight units arrive in a given area with a probability P at least 
80%. 

Solution. According to the graph Fg(x) aaa P = 0.8, the quantile is found X9g = 3.1. 

Since x = (t —m;,)o™}, then t = m, + 0,X9g = 3 + 0.5 * 3.1 = 4.55. 

It is necessary to dwell on one case of asymptotic behavior of the hypernormal distribution. 
Let n > © (practically for x > 3), a natural consequence of this condition is F,(x) > 1. Then 
the differential equation (1) can be represented for standard conditions (m = 0,0 = 1) in the 
following form: 

nF, (x) + xF,(x) = 0. (33) 


Separating the variables, it is found: 


dln fn(x x 
dln fn(x) Sy (34) 
ax n 
__ aFn(x) 
TAe f n(x) = are ae — MAOTHOCTS IMMepHOpMaAbHOTO paciipeAcAcHHA. 


Integration of equation (34) makes it possible to verify the validity of the following result, 
presented in the form of the following theory. 

Theory 2. For large values of the argument, the hypernormal distribution asymptotically 
tends to the normal distribution with density (5). 

The result obtained can be somewhat strengthened by considering, instead of the 
hypernormal distribution function, the normal distribution function that satisfies the differential 


equation (1), with a variance depending on the value of the argument: 
2 


x 
1 De. 


fn(x) = (2D (x) 


Using equation (1), it can be shown that the nature of the change in the variance D(x) is 


determined by the following differential equation 
(x? -D)=- 2xD ee 0. (35) 
n[p(x)] 7 
where D(x) is Laplace function. 
It follows from differential equation (35) that for 
x > 0,D(x) > 02 n. 


We set ourselves the goal of determining the function of the limiting hypernormal 
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distribution (extreme distribution of type 2) in the following form: 
F(x) = P{VnXx < x}. (36) 
For sufficiently large values of n. For standard conditions (m= 0,0 =1), the 
differentiated equation (1) can be represented (apu n — 0) in the following form: 
NF (X) Foo (X) + XFio(x) = 0. (37) 
It is possible to check that the change of the independent variable y = = allows to 
transform equation (37) to the following form: 


2 
2F.,(y) a2 Foy, Feo 


= — = 0.135) 
dy dy 
Separating the variables and integrating, it is found: 


ee = | 9) 
dy 2 
From this it follows: 
P dFoo 
y= —1 0 In Foo = —2E; (In P) (40) 


where 1;P and E; (InP) is integral logarithm and integral function respectively. 
Making the reverse transition from y to X, it can be obtained the following result, presented 
in a compact form in the form of the following theory. 


Theory 3. The quantile function of the hypernormal distribution asymptotically (n > 0) 


Xp = V2n/—E,(nP). 
The proof of the theory follows from the above reasoning and confirms the validity of 
relation (6). 


approaches the function 


As an example, illustrating the applicability of an extreme type 2 distribution, consider a test 
planning problem. 

Example 3. Tests for the failure-free operation of a product should be performed for no 
more than 100 days. It is assumed that by the end of the tests at least 90% of the ordered 
products should fail. A preliminary assessment of the reliability indicators showed that the 
operating time is average, but the failure T, is 50 days, and its standard deviation 0, is equal to 
5 days. How many items do you need to order and test? 

Decision. If the law of the time-to-failure distribution is not known, and the ratio between 
the mathematical expectation and the standard expectation is such that there are grounds to 
consider the range of permissible values practically unlimited and the estimated number of 
ordered products n > 10, then using the function of the limiting hypernormal distribution, 
according to Table 2, we define the argument (Table 2): 

F,.(%) = 0.9 (% = 1.9). 


Then: 
Tg7To _ 
he 1.9 
where Tg = 100 days (directive test time). 
Therefore: 


Tg-To \2 =BO\ . 
n= (=) = (=) =~ 25 items. 


071.9 5+*1.9 
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When solving practical problems, one can hardly hope for the completeness of information 
about the initial random variable, which makes it possible to estimate higher moments of its 
distribution. If the available information about the initial random variable allows us to give only 
an estimate of its mathematical expectation, then in the conditions of problem (24)-(29) the 
condition (28) is excluded, as will be shown below, the entropy maximum is achieved on a 
distribution that satisfies the differential equation (7). By inverting the distribution function 
G, (x), one can obtain the following result, presented as a theory. 

Theory 4. Let X be a random variable with density g(x) > 0,x € (0,0), G(x) is the 
distribution function of the maximum value from the set of n random variables from the general 
population defined by the random variable X. Let, further, the mathematical expectation of the 
random variable X: 

m= f, xg(x)dx. 

Then the entropy maximum is reached on a distribution that satisfies the differential 
equation (7) with the quantile function (9). 

The proof will be as follows. Peers equation (7) is the Euler-Lagrange equation: 

n[Gy (x)] = mes (x) — vG, (x) = 0 (41) 

G,(0) = 0,G,(c) 

where v is undefined multiplier. 
The next variational problem: 

=— a G,, (x) In G, (x)dx > max, 

te G,(x)dx = 1, 

J, g@)dx = 1, 

J, xg(x)dx =m 


Gé.G@= LU gGddx]". (42) 
Formal integration of equation (41) for n = 1 gives an exponential distribution law, and for 
n # 1 leads to a dependence of the form 
x = fo" — Aon (43) 
vG4Gq(0) 
1 


Using the substitution Z = Gr (43), upeoOpasyem kK TAaOATHHOMy: 
L 
i ACL Ba) Van = = 2F\(v,u; 1+ uw; —Bu) 


where 2F,(a,p;y;S) is hypergeometric Gaussian function. Therefore: 
Gy (0)x = Z” »F,(1,n;n + 1; —Bz) (44) 


ot 
G,,(0)x = poF,(1, nyn+ 1;B%/p), 
where B = 5 a up=G,(x)= Pix (") x}. 
The normalisation condition for the distribution function G,(%) implies the following 
relation: 
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als 
Indeed, the quantile functions ne from relation (44) have the form: 


=5 aT 2F,C1,n;n + 1;-BZ). 


For the original random variable: 


Xp» = = 2F,(1,n;n +15 —B"/p). 


n . 
For the maximum value from a set of n random variables. For X, > 00,Z = Gn (x) > 


1 and p = G,(%) > 1 and vice versa. 
The Gauss hypergeometric function can be represented as the following series: 


co a4rBrs™ s 
iGO Vis) A205 
Yor) 7 
wCaes 
I'(a) 


For p = 1(z = 1), series 


where A(,) = 


o lmnr@ (- oe 
ae a 


should diverge. 
According to the d’Alembert test, series (45) will diverge if the relation follows: 
LereayMrsiy(—BYOTVn+D (yr! I 
4+n)r+aylaynay(-B)"(r+1)! 

Or after transformation: 

_ B = ntrt2 > 1. 
“n+2— 

Therefore, in order for the series to diverge (the normalisation condition is satisfied), it is 


necessary that for the sufficiently large r(r > ©) complex: 


On 
To determine the density value at the initial point G,,(0), one integrates relation (44) taking 


into account the obtained value for {8: 


e G,,(0)xdz = hyn F,(1,n;n + 1; z)dz. 
0 0 2 


The left side of the equation in accordance with the definition z = [Gy (xyz is the product 
of the value of the distribution density at zero G,,(0) and the mathematical expectation of the 
original random variable. It can be shown that the right side of the equation is equal to one. 
Indeed, by representing the hypergeometric function as a series and changing the order of 
summation and integration, it is defined: 


i En TOMO" gy = ye LOG) Fa n 
r=0 (+n), 1! r=0 (440); 1! (n+r)(nt+r+1)’ 


Series })7-9 ——— can be represented as the difference between two series: 
nas 


n 


co 1 
Dro Ganeeney = PLre0 Gary ~ MLro0 Gap 


which after transformation can be represented as follows: 
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1 1 1 d 
Dr 6s on nYir=0 —1NYm=nt1 ae dim=n = 


aarp 


Thus, G,(0)m = 1. Therefore, G,, (0) = 


=a: 


and the Lagrange multiplier v = 6 G ~(0) = 
—m 

Substituting the value of the Lagrange multiplier into the differential equation (41) makes it 
possible to verify the validity of the differential equation (5) and the quantile function (9). The 
theory has been proven. 

It is also necessary to clarify the method for determining the mathematical expectation of 
the largest value for the case under consideration. In a way analogous to that which was applied 
in the derivation of the formula for determining the value of the density at zero G,(0), one can 
show that: 


| = 00 larry) 
E[Xy = dr=0 (n+r)(2n+2)" 


To sum up, it suffices to integrate the quantile function (9): 


1 1 3 
E|x| = J, xpdp =m J, pF, (1,n;n + 1;%/p)dp. 


Representing the hypergeometric function as a series and rearranging the operations of 


summation and integration, one can find: 
2 


E|x| =m 6 one fp Lndp = mY2. TCEES OT EoT 

In table 3 shows the values of the sum S,, which make it possible to estimate the 
mathematical expectation of the maximum value from samples of n to 10 (for n > 10, one can 
use, as will be shown below, the asymptotic properties of the obtained extreme distributions of 
extreme values) (Table 3). 

If we set ourselves the goal of determining the function of the limiting extremal distribution 
of type 4, then, by performing constructions similar to the constructions used in the proof of 
Theory 3, we can obtain the result presented in the form of the following theory, given by virtue 
of obviousness without proof. 

Theory 5. The quantile function of type 3 extremal distribution asymptotically (n > 0%) 
approaches the function: 


x = —7 FE, (InP). 
In—— 


As an example, illustrating the applicability of extremal distributions of types 3 and 4, 
consider the problem of “forecast by one point”. 

Example 4. A discrete random process X(N) is observed, its value x(1) is fixed at the first 
observation point. What is the expected value of the maximum value at the second point at the 
20th pointe 

Solution. It is natural to take m = x(1) as an estimate of the mathematical expectation of 
the average process under one observation. Then, using Table 3, one finds E [max X (2)| = 


x(1) * 1.1628 and, using dependence (14), one determines: 


E [max xX(20)| = x(1) id = x(1) * 14.2067. 
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Further, it is necessary to consider some features of the construction of extremal 
distributions of minimal values. Differential distribution type 5 Q,(%) is the Euler-Lagrange 
equation of a variational problem similar to problem (24)-(29). The difference in the formulation 
of the variational problem lies in the replacement of relation (29) by the dependence: 

On(x) = 1-[1— FF" (46) 
defining the distribution function of the smallest value Q,(%) through the distribution function 
of the original random variable F (x). This implies that 
He = —JS> Qn(x) InQ,(x)dx > max 
under isoperimetric conditions (25)-(28) and holonomic constraint (46). 

The solution of this variational problem allows us to formulate a result similar to that 
stated in Theory 1. 

Theory 6. Let X be a random variable with density f(x) > 0,x € (—0, 00), Q,, (x) is the 
distribution function of the smallest (leftmost) member of the variational series constructed from 
a sample of finite size n from the general set determined by the random variable X. Let, further, 
the two first central moments of the random variable X: 

m= |" afaax, 
o2 = [-@ —m)*f (x)dx. 

Then the entropy maximum is reached on a distribution that satisfies the differential 
equation (15). 

Using dependence (46) and differential equation (15), one can verify that the extremal 
distribution Q,(x) corresponds to the distribution function F(x) of the original random 
variable, determined by differential equation (3). 

As an example of using a type 5 distribution, consider the following problem. 

Example 5. Under the conditions of Example 1, find the probability of preparing the first 
product by the time T = 100 min. 


Solution. The solution to this problem is reduced to a sequence of reasoning and actions 


applied in solving Example 1. As a result, passing to the normalized standard values x = —= = 
= = 0 according to the schedule Q;(x) for x = 0, one finds Qs (x) = P(tay < 100) = 
0.72. 


Analytic properties of the function Q,(x) are similar to those of the function F,(X). 
Therefore, for small values of the argument x <m-— 30 the extreme distribution of type V 
asymptotically approaches the normal distribution about the parameters m and no, and for 


n — 00 it degenerates into an extreme distribution of type 6, the quantile function of which is 
described by dependence (17). 


Generalising the results concerning extremal distributions of types 1 and 5, it seems 


appropriate to find the distribution function of the order statistics me ) (m = 1,2,...,n), which 


provides the entropy maximum under isoperimetric conditions (26)-(28). If the initial population 


distribution F(x) has density f(x), then the distribution of order statistics x ) umeer 


TIAOTHOCTb BY Aa: 
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fam (20) = = P{XK” < x} = F(xc)™-*[1 — FOX) -™ F(x) (47) 


The seen distribution of the order statistics (the favorite member of the variational 


(m— rf ee m) 


series) is determined by the solution of the Euler-Lagrange equation of the following variational 


problem: 
=~ Jy fn) In frm (x) dx > max (48) 
1 = fro, fam (x) dx (49) 
m = J", xf (x) dx (50) 
o? = f° —m)? f(x)dx (51) 


! F(x) 
Fum(x) = ena fe iim @ Rea) amar va (cyA 


Applying the course of reasoning and those transformations and constructions that were 
used in the proof of Theory 1, we can obtain a nonlinear differential equation with respect to 
the function Fim (x): 

z (x-m) D225 [Fnm (x)-F(xo)l” DY’Fnm 
Z vl 
OL 
ibe Fnm (Xo) ( ) 
Satisfying the boundary conditions Fy (—%) = 0, Frm(co) = 1. 


_ (m- mine m)! 


Depending on (53) D FI™1-F ynon = is conversion operator with 
properties DY = DCDY-*),p?'= B, 

The reference value X9 for the function Fy,(x) can be chosen at any point where the 
distribution density (47) does not vanish. Note that for m= n differential equation (53) 
degenerates into equation (1), and for m = n into equation (15). 

Thus, the solution of the variational problem (48)-(53) and the representation of the 
nonlinear part of the differential equation using an operator series (S. Lie’s series) allows us to 
represent the result obtained in a compact form in the form of the following theory. 


Theory 7. Let X be a random variable with density f(x) > 0,x € (—©, 00), Fim (X) is the 


‘s ‘ 7 a a ‘é n ‘ * = 
distribution function of the order statistics x”) for a variational series constructed from a 
sample of finite size n from the general population determined by the random variable X. 


Further, let the first two central moments of the random variable X be known: 
m= (jee xf (x)dx 
o* = (ie —m)*f (x)dx. 


Then the entropy maximum is reached for a distribution that satisfies the differential 
equation (53). 

In conclusion, it should note some features of extremal distributions of types 7 and 8. 
Differential equation (18) is the Euler-Lagrange equation of a variational problem similar to 


problem (42). The difference in the formulation of the problems lies in the replacement of the 
last holonomic constraint G,(x) = [ f, g (x)dx]" by the dependence R,(x) = 1- [1 _ 
i r(x) dx]", which determines the distribution function of the smallest value R,(x) through 


the distribution function of the original random variableR(x) = f ° r(x)dx. 
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Integration of the differential equation (18) allows one to obtain the quantile function (20), 
and the passage to the limit n > 00 yields dependence (22). A characteristic feature of these 
distributions is that they are realised in the class of finite functions (truncated on the right). It 
should also be noted that the dependencies for the mathematical expectations of the smallest 
values (21) and (23) can find practical application in the express evaluation of samples from an 
unknown general population. 

The asymptotic behavior of the largest observation in a sample of size n from a distribution 
with distribution function F(x) was a problem in the “classical” theory of extreme values. The 
central result of this theory (the theorem on three types of limit distributions) was first obtained 
by Fisher and L. Tippett in 1928 (T7ppezt, 2072) and was later proved in full generality in 1943 by 
B.V. Gnedenko (Gnedenko, 1943). A systematic exposition of the theory of limiting distributions 
of extreme quantities and applications to technical problems can be found in Gumbel’s 
monograph (Gumbel, 2012). However, here it is necessary to note some specific limitations of 
the asymptotic theory of extremal quantities. First of all, all extreme value distributions are 
derived either from the exact original distribution or from a distribution of some type. The 
original distribution from which the extreme values are selected must belong to one of the three 
types of distributions. In reality, the analytic properties of the original distribution are rarely 
known, and hence the conditions for using the asymptotic theory of extreme values do not 
always correspond to observations and practical applications. Note that in the “classical” theory 
of extreme values, when constructing parametric forms of distributions of extreme values, the 
key idea (stability postulates) used earlier by Fisher and L. Tippett is used, which consists in the 
following. Since the largest observation in a sample of size mn can be considered as the largest 
member in a sample of size N, consisting of the maximum members of samples of size m, and 
since in the case of the existence of a limit distribution A(x) both of these distributions will tend 
to A(x) at m > 00, then A(x) must satisfy the relation A" (a,x + b) = A(x), ie., the largest 
observation in the sample of size n from the distribution with the distribution function A(x) 
must, after appropriate normalisation, itself have a limit distribution function A(x). The solution 
of this functional equation with respect to A(x) allows us to obtain the following three 
parametric forms (three types of distributions of extreme values built on the stability postulate): 

Type 1: Ay (x) = exp(—e~*) , -—0 < xo, 

0,x <0 
Type 2g) = = (-x7*),a>0,x >0° 


exp(—(—x)%),a>0,x <0 
Type 3:Ag(x) =f p(-( eae . 


In this case, it is logical to compare the distributions of extreme values based on different 


construction principles. As an example, consider the graphs of distribution functions F,,(x) and 
A, (x) Ay(x) = exp (-e > ). Naturally, the functions must be compared with the same scale 


and position parameters. Note that if a sample X1,%2,...,X_ is given from the population with 
= 2 . _ x-a . 

the distribution function A,(x) = exp (—exp (— aE as estimates a* and b* of unknown 

parameters a and b, one can take solutions of equations (6): 


lok x;-a* = 
pi=t b* =C 
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lor xj-a*) 
k faf{- b* j=, 


where C is Euler constant. 


= 1 
Equating the empirical values of the moments of random variables pit x; and 


pyr exp (-=) with the theoretical ones for the distribution F(x), after some algebraic 


transformations, one can find equations relating the parameters a* and b* of the distribution 
function A, (x) and m, and oV2n of the distribution function F(x): 
m+oav2nl, = a* + cb* (54) 


== ov2n ee 
e bl, =e b 
b 


whete 


I, = J) —E;(np)dp = 0.6967, 
2 1 _aV2n,/—EFi(in p) 
EG ()=s e P dp. 


With appropriately chosen normalizing and centering parameters a* = 0 and b* = 1, the 


limiting hypernormal distribution corresponds from equations (54) to the values of the 
parameters oV2n and m, determined from the following relations obtained from the systems 
of equations (54): 
oV2n +1n I, (oV2n) = 0, 
m =0-oVv2nI,. 

The figure (Figure 21) shows the graph of the function A(x) = e-° ” and the graph of the 
function of the limiting hypernormal distribution corresponding to it in terms of parameters 
F(x) = -1.48 + 3,/—-E,(In P). 

A comparison of the graphs indicates the closeness of the statistical laws of extreme values 
obtained under different assumptions. On the other hand, it is not possible to introduce partial 
ordering into the set of distribution functions of the form A, (x) and F,,(%) (the graphs illustrate 
this circumstance): the function A,(x) does not dominate the function in the sense of the first 
otder. Therefore, when choosing a mathematical model of the mechanism for the formation of 
extreme random variables, when constructing estimates that guarantee their effectiveness, etc., 
one should be guided by the stability postulate, if it is known that the distribution of the initial 
random variable belongs to the exponential type distribution (4), and by the principle of 
maximum uncertainty, if only estimates of the mathematical expectation and the variance of the 
original random variable (1) are known. 

As an example, illustrating the possibility of using the considered laws of distributions of 
extreme values, it is presented the distribution function of the levels of water rise at the mouth 


of the Neva, constructed according to statistical data fixed since 1703 (7): 
X-221,4 
A,(x) = exp {=e 30,38 \ 
x = 183 + 7.97V2n,/—E, (In P) 


where x is level of Neva’s mouth water rise in centimeters. 
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Comparison of the probabilities of occurrence of events P{x < 200}", built on these 
dependencies, is similar to the nature of dependencies (Fzgure 27). 

The discussion of the applications of the theory of extrema of sequences of random 
variables is given fragmentarily and does not yet cover the whole variety of possible problems 
related to extreme values, because “... nature speaks to us in the language of mathematics” (G. 
Galileo) and it is quite appropriate that this area is worthy of study on its own. yourself. 
Therefore, it seems appropriate to proceed to the consideration of the basic problems of 
constructing models and some generalisations, confining ourselves mainly to the formulation of 


the main results. 


Conclusion 

Thus, in the study course, the definition of an extreme value within the framework of the 
theory was refined, the typology of the distribution of maximum values was analysed, seven 
theories of the hypernormal distribution were identified and their proofs were presented, and 
practical examples of the application of each theory were given. The practical significance of the 
study of extreme random variables models in various areas of industrial human activity was 
confirmed. 

The materials of the study can be used in the widest range: from application in risk 
management of industrial production to predicting the probabilities of natural phenomena, 
which makes it possible to prevent significant economic and social losses of society, as well as 
make a tangible contribution to programming the probabilities of the development of the society 
of the future. 


Do 
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Appendix 


Table 1. Extreme distributions and statistical characteristics of the original random variable 


Determinatio Extreme distributions and statistical 
ata characteristics of the original random variable 
Matieananioal Sxpeciion m and Expected value m 
variance 0 
PEX®) < x} F(X) F(X) Gy (X) Goo (X) 
type | type 2 type 3 type 4 
(limiting (n > 00)) (limiting (n > ©)) 
Qn(X) Qo (X) Ry (X) R(X) 
type 5 type 6 type 7 type 8 
P{X’ > x} (limiting (n > ©)) (limiting (n > )) 


Table 2. Values of the function of the limiting hypernormal distribution F(x), whose argument is the 


value X 

x F(X) x Foo(X) 
0,1 0,0230 1,6 0,830 
0,2 0,0707 1,7 0,858 
0,3 0,126 1,8 0,873 
0,4 0,192 1,9 0,904 
05 0,259 2,0 0,921 
0,6 0,326 a4 0,936 
0,7 0,393 22 0,949 
0,8 0,457 5.5 0,959 
0,9 0,518 24 0,967 
1,0 0,573 25 0,975 
11 0,628 26 0,981 
12 0,676 27 0,985 
1,3 0,722 28 0,988 
1,4 0,762 29 0,992 
15 0,798 3,0 0,996 


Table 3. The value of the sum S,, of a series from 1 to 5 


n 1 2 3 4 5 

Sn 1 1,628 2,2643 2,8878 3,4967 
n 6 7 8 9 10 
Sn 4,0904 4,6694 5,2338 5,7844 6,3214 
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Figure 1. Distribution function F,(x) = 
Py xt < x}; (m = 0,07 = 1) 


Figure 2. Distribution function F,(x) = P{x? < 
x}; (m = 0.4634; a? = 1.1077) 
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Figure 3. Distribution function F3(x) = 
P{x® < x}; m = 0.5355; 0? = 1.3594 
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Figure 5. Distribution function F,; (x) = 
P{x® < x}; m = 1.1355; 0? = 2.0190 
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Figure 4. Distribution function F,(x) = 
P{x® < x};m = 0.9764; a? = 1.6622 


x 
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Figure 6. Distribution function Fy (x) = P{ x < 


x}; m = 1.2458; o? = 2.3316 
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Figure 7. Distribution function F;(x) = 
P{x < x}; m = 1.4656; 0? = 2.9520 
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Figure 8. Distribution function Fg(x) = 
P{x® < x}; m = 15504; 0? = 3.3127 
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Figure 9. Distribution function Fy(x) = 
P{x® <x}; m = 1.5748; 0? = 3.5122 
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Figure 11. Distribution function Q,(x) = 
PL < x}; m=0;07=1 
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Figure 10. Distribution function Fi (x) = 
ccm < x}; m = 1.5792; 0? = 3.6438 
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Figure 12. Distribution function Q)(x) = 
P{x® < x}; m = -0.5321; 0? = 1.0492 
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Figure 13. Distribution function Q3(x) = Figure 14. Distribution function Q,(x) = 
P{ x6 ig x}; m = —0.7330; o2 = 1.1250 P{x6 4) < x}; m = —0.9445; o2 = 1.1546 
Qs 
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Figure 15. Distribution function Q,(x) = Figure 16. Distribution function Q(x) = 
P{ x6 (5) < x}; m = —1.0692; 02 = 2.1477 Pix! (6) ~ x}: m = —1.1950; o2 = 2.5258 


Figure 17. Distribution function Q7(x) = Figure 18. Distribution function Qg(x) = 
P{ x6  < x}; m = —1.2977; o2 = 2.9194 P{ x6 (@) < x}; m = —1.4451; 02 = 3.2427 
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Figure 19. Distribution function Qg(x) = Figure 20. Distribution function Q)9(x) = 
P{x® <x}; m = -1.5573; 0? = 3.6019 P{x0 < x}; m = 1.5703; 0? = 4.0949 
F(x) 
1,0 
0,5 
=2. -148 -1 0 1 2 3 
—-—-—— - A, (x)= 07%" 
— XY, =-148+3 8, (InF) 


Figure 21. Distribution function F,, (x) * A, (x); m = 
—1.5573; 0? = 3.6019 
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